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Abstract 

^ ' The validity of the approximations done in the mean field description of 

P3 . the quasi-elastic excitation of medium-heavy nuclei is discussed. A test of 

the reliability of the plane wave Born approximation is presented. The un- 
certainty related to the choice of the electromagnetic nucleon form factors is 
discussed. The effects produced by the meson exchange currents generated 
P^ \ by the exchange of a single pion are studied, as well as the need of including 

relativistic corrections to the one-body currents. The results of the continuum 
shell model and those of the Fermi gas model are compared and the need of 
treating the emitted nucleon within a relativistic framework is studied. We 
analyze the results of Random Phase Approximation calculations not only in 
terms of their effects on the mean field responses, but also in terms of the 
theoretical consistency in the choice of the effective residual interaction. The 
role of the final state interaction is investigated. A microscopic justification 
of the use of complex optical potential to describe the final state interaction 
is provided. The effects of the final state interaction on the mean field re- 
sponses are studied with a sum rule conserving model. A comparison with 
the experimental data measured in the ^^C and ^''Ca nuclei is shown. 
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1 Introduction 

The Mean Field (MF) model is the basis of any description of many-particle sys- 
tems. In this model the particles composing the system are supposed to move 
independently from each other in an average potential. Under this assumption the 
complicated many-body problem is transformed into a set of, easy to solve, one-body 
problems. The solution of the MF equations produces a set of orthonormal single 
particle (sp) states whose tensor product forms a basis in the many-body Hilbert 
space. This basis is used in more complicated many-body treatments. 

Usually the behavior of the real many-particle systems is quite different from the 
MF predictions. There are however phenomena which can be rather well described 
by the MF model. In these cases the more complicated many-body effects can be 
treated as a correction to the MF solution. One of these phenomena is the quasi- 
elastic (QE) excitation which is dominated by the single particle dynamics, well 
described by the MF approach in terms of one-particle one-hole (Ip-lh) excitations. 

Experimentally, the QE excitation is obtained by using probes which interact 
weakly with the system and by controlling the kinematical conditions in such a way 
that only one of the particles composing the system is directly struck by the external 
probe. This produces a relatively small perturbation of the system which however 
responds as a whole because the knocked particle is linked to the other ones. The 
total response is related to the characteristics of the interaction among the particles. 

In the electron excitation of the atomic nuclei, the region where the nuclear 
response is quasi-elastic is well identified. The inelastic electron scattering cross 
section shows an evident peak at higher energies than those of the giant resonance 
region. The position of this peak is located at values very close to |qp/2M, where q 
is the momentum transferred from the electron to the nucleus and M is the nucleon 
mass. The width of the peak increases when the mass number of the target nucleus 
increases. These two facts lead to the interpretation of the peak as produced by the 
elastic scattering of the electron with a single nucleon. An ideal experiment on a free 
nucleon produces a narrow peak positioned exactly at |qp/2M. Understanding the 
deviation from this situation means understanding the nuclear electron excitation 
mechanism and the nuclear many-body dynamics. 

The first electron scattering experiment aimed to probe the nuclear QE response 
was done at the Stanford Mark III linear accelerator in the early seventies ||Mon71|| 



The data were analyzed within a Fermi gas (FG) model ||Mon69|| , i.e. considering 
the nucleus as an infinite system of non interaction fermions, and the corresponding 
fits were obtained by changing the two parameters of the model: the average nucleon 
interaction energy, {E), and the Fermi momentum, k-p. The position of the peak is 
fixed by {E), while both height and width are determined by the Fermi momentum. 
The interesting result was that the three quantities could be reproduced so well by 
two parameters. 

The aim of this experiment work was to obtain an empirical estimate of the value 
of k-p. The commonly accepted value of kp (about 270 MeV/c) is obtained from the 
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Figure 1: FG model results compared with the experimental data of refs. [Yat9S, Wil97] 
for the ^''Ca nucleus. The electron energies e and the scattering angles 9 are indicated in the 
figures. The dashed lines show the results obtained by using the values of /cf=251 MeV/c 
and {E)=28 MeV fixed in ref. | Mon71 |. The full lines have been calculated changing 
the values of the Fermi momentum and of the average binding energy to reproduce at 
best the experimental data: (a) kp=340 MeV/c and {E)=0 MeV, (b) /cf=330 MeV/c 
and {E)=5 MeV, (c) /i:F=330 MeV/c and {E)=5 MeV and (d) A;f=320 MeV/c and 
(^)=90 MeV. 



empirical value of nuclear matter density, p = 0.17 fm~'^, fixed by considering the 
nucleus a homogeneous sphere of radius R = tqAs, with tq = 1.18 fm given by 
elastic electron scattering experiments. The fit procedure used to reproduce the QE 
data provided for heavy nuclei a roughly constant value of kp of about 260 MeV/c, 
in good agreement with the value obtained with the usual procedure. 

The success of this description of the data within a FG model suppressed inter- 
est in the investigation of the QE peak for about ten years, until the data of the 
separated charge and current responses became available ||Alt80|| . The FG model 
which successfully explained the cross section was unable to describe separately its 
two components ||Alt8[ 



Looking back at the analysis of the pioneering Stanford experiment we find that 
what is nowadays considered interesting physics is hidden in the two parameters of 



the fit. The nuclear structure problem and the interplay between single-particle and 
many-body dynamics have been obscured by the fit procedure. 

In fig. |I| we compare the FG results with the recent ^°Ca data taken at Bates 
yat93| , [Wil97|| . The values of {E) and k-p defined in ref. |[Mon71|| produce cross 



sections well above the data (see dashed curves). On the other hand, if {E) and k-p 
are fixed to reproduce one data set (full lines) the other ones are not well described. 

The FG model provides a qualitative description of the phenomenon, but to 
obtain a quantitative description of the different sets of data it is necessary to go 
beyond this model and treat explicitly those degrees of freedom that are hidden in 
the effective values of the two parameters of the model. 

The aim of this article is to discuss the successes and failures of the MF de- 
scription of the electron induced QE excitation of medium and heavy nuclei. The 
discussion will be focused on what can be learnt about the interaction between the 
electromagnetic probe and the nucleus and the dynamics of the nuclear excitation. 
In the presentation we shall try to separate the problems concerning the reaction 
mechanism from those regarding the nuclear structure. This separation is somewhat 
artificial but it is useful to better determine and clarify the various approximations 
and uncertainties done in the study of the QE peak. 

The description of the reaction mechanism is based upon the Plane Wave Born 
Approximation and the use of the traditional electromagnetic one-body operators, 
and we shall discuss the limits of all these approximations considering the distortion 
of the electron wave functions, the presence of Meson Exchange Currents and the 
relativistic corrections to the one-body operators. 

The nuclear structure problem is tackled within the non-relativistic formalism. 
The basic model we shall adopt is the Continuum Shell Model, where the nucleons 
move independently from each other in a finite size average potential. We shall 
present, however, results obtained by describing the nucleus as an infinite system of 
non interacting nucleons. In addition to these MF descriptions of the excitation we 
shall also discuss results obtained with the Random Phase Approximation theory. 



2 The electron-nucleus interaction 

Our understanding of the electromagnetic processes is certainly deeper than that 
of those processes dominated by the strong interaction. There are however various 
problems not very well solved, or defined, in the description of the interaction be- 
tween electrons and nuclei. In this section we shall discuss some of them. We shall 
present first what we call the standard derivation of the inelastic electron nucleus 
cross section. Then we shall discuss the validity of the approximations done in this 
approach, with special attention to their influence on the study of the QE excitation. 
In particular, the role of the nucleon form factors, the relativistic corrections and 
the Meson Exchange Currents are discussed. 



2.1 The standard treatment 



In this section we present the basic points of the usual derivation of the electron 
scattering cross section off atomic nuclei. More detailed presentations of this deriva- 
tion can be found in various review articles (see, e.g., ||deF66| , Pon75| , |Cio8CI| , Pre89| , 
BoT5§). 

Throughout the paper we shall use the Bjorken and Drell ||Bio64|| metric and 
the convention that the repeated indexes of a four-vector imply a sum on them: 
afj,b'^ = aolf — a ■ b. 

Using the standard Feynman rules, the scattering amplitude for a single photon 



exchange process can be written as ||Bio64|| : 

{f,kf\S\i,h) = -t fd'x[ei^f{x)-f^A^{x)ij,{x)]. 



In this equation the initial (final) state of the system is defined as the tensor product 
between the initial (final) state of the electron and of the nucleus. As we work in the 
Plane Wave Born Approximation (PWBA), the electron wave functions are plane 
wave solution of the free Dirac equation: 



m 



^{x) = J^u{k,s)e-"'"''^ 



(2) 



where u is one of the positive energy four-component Dirac spinors | Bjo64| |. The 
electron states are characterized by the values of the four-momentum k = {e, k}. 
The final formulae will be independent of the normalization volume V. 

In eq. (|lD we have indicated with A'^ the classical electromagnetic field generated 
by the nucleus. This should satisfy the field equation: 



OA^ix) = e(/|J^(x)|2) 



(3) 



where D = V^ — d'^/dt^ and J^ is the nuclear current operator. The calculation now 
proceed by inserting the solution of eq. (^ and the expression of the electron wave 
functions (|^) in eq. (|I]) which can be rewritten as: 



{f,kf\S\i,ki) = -ie^{kf\j^\ki) 



-^fii/ 



q^Qp 



where the four-momentum transfer is defined as g^ 
current is given by: 



d'xe-'^"""(/|J,(x)|i) 



(4) 



fcj and the electron 
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(5) 



The nuclear current is supposed to be described by an operator local in r. Its time 
dependence is given by: 



Mt,r) 



JHt 



J^^ir)e 



-iHt 



(6) 



After performing the integration on the time variable the scattering amphtude can 
be expressed as: 

{f,kf\S\t,h) = -^e'27T5{E,-Ef-Lu){kf\J,\h)^{f\J'^{cl)\^), (7) 

where the Ei and Ef indicate, respectively, the initial and final energy of the nucleus 
and uj is the energy transferred by the electron. The expression (0) depends explicitly 
on the Fourier transform of the nuclear current 

J^(q) = Jd'xe^'^-^J.ir). (8) 

The nuclear current should satisfy the continuity equation: 

9^J^(x) = 0, (9) 

which, using the exphcit time dependent expression (^, can be written, in coordinate 
space, as : 

V • J(r) = -iiuj\r) (10) 

while in momentum space it assumes the form: 

q-J(q) = ujp{q), (11) 

where we have used J° = p. The continuity equation imposes a relation between the 
four components of the nuclear current, therefore only three of them are independent. 
It is necessary to choose the independent components to be used to express the 
scattering amplitude. It is convenient to work in a reference system where the 
quantization z-axis points along the direction of the momentum transfer q. Eq. ([TT|) 
shows that the component of the current along this axis, the longitudinal component 
Jz, can be expressed in terms of the charge operator p. Therefore the scattering 
amplitude can be expressed in terms of p and of two components orthogonals to 
the direction of q. These are called transverse components and they are usually 
expressed in spherical coordinates: 
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T^{J,±i3y). (12) 



To calculate the cross section one has to square the scattering amplitude, to 
average on the initial states and to sum over the final ones. Since the electron 
energies e are much bigger than its rest mass m, it is common to use the so called 
ultrarelativistic approximation consisting in neglecting the terms depending on m. 
The calculation is rather long but it can be done without any further approximation. 
The cross section assumes the expression: 



d^ = "Mo.. I (^)'fl.fe,-) + (tan= 5-^1 flr(9,-) } , (13) 



^)flH.,-) 



where 9 is the scattering angle, defined as the angle between the directions of q and 
k/, 9 = |q| and UMott is the Mott cross section 

_ / «cos(g/2) V 
where we have indicated with a the fine structure constant. 



In the expression ([T3| ) the cross section is composed of two parts, each of them 
formed by a kinematical term and a response function containing the full nuclear 
structure information. The longitudinal and transverse response functions are re- 
spectively defined as: 

RL{q,uj) = Y.T.^(Ef-E,-uj)\{f\piq)\i)\' (15) 

i f 

Mq,^) = T.ESiEj-E.,-co){\{f\.MqM' + \{f\J4qM'). (16) 

i / 

The bar on the sum on i indicates the average on the initial states. 

If the expression (O) of the cross section is valid it is possible to obtain an 



experimental separation of the longitudinal and transverse responses. To achieve 
this goal it is necessary to perform a set of measurements at fixed u and q but for 
different values of the scattering angle 9. The ratio of the cross section and UMott 
should be plotted against vt = tan^(6'/2) — (tu^ — q'^)/(2g^) and all the points should 
lie on a straight line. The slope of this line gives Rt while the intersection with the 



y axis provides Rl- This procedure is called the Rosenbluth separation [[Ros50 



The nuclear responses are calculated under the hypothesis that the nucleus makes 
a transition from one state of definite angular momentum \i) = \JiMi) to another 
one I/) = \JfMf). For this reason it is convenient to rewrite eq. (^ by making 
a multipole expansion of the exponential. Considering the parity of the nuclear 
final state, it is possible to separate the electric and magnetic transitions and at the 
end one has to deal with three type of multipole operators: the Coulomb operators 
producing transitions of natural parity. 



Mjuiq) = J d'rjj{qr) YjMif) p(r) , (17) 

the electric current operators which give rise also to transitions of natural parity, 
TfMiq) = -Jd'r{Vx [jAqr)Yfj{f)]} . J(r) , (18) 

q '' 

and the magnetic current operators producing transitions of unnatural parity, 

Tfuiq) = Jd'rjj{qr)Yyj{n).3{r). (19) 

In the eqs. (|T^-(|TU|) jj indicates a spherical Bessel function, Yjm is a spherical 
harmonics and Yj£ is a vector spherical harmonics ||Edm57| . 



We shall deal with scattering from doubly closed shell nuclei, therefore Jj = 
and then the value of the angular momentum of the final state correspond to the 
multipolarity of the transition. The nuclear responses can be written as: 

oo 

RL{q^u) = 47r^^ |(/5,J||Mj(g) ||0)|2 5(E^-^), (20) 

/3 J=0 

oo 

RT{q.u) = 4vrE5: {\{f3, J \\ Tf{q) \\ 0)\'+ 



\{f5,J\\Tf{q)\\Q)\')5{E,-u), (21) 



where 13 indicates all the quantum numbers characterizing the nuclear final state 
other than J. 

The only elements which remain to be specified are the explicit expressions of 
the nuclear electromagnetic operators p and J. Up to now the only hypothesis made 
on these operators has been done in eq. (^. Since in general the nuclear many-body 
wave functions are described using nucleonic degrees of freedom in a non relativistic 
framework, the standard treatment considers the electromagnetic operators pro- 
duced by non relativistic pointlike nucleons. This means that the nuclear charge is 
composed by the charge of the nucleons and the current is composed by the move- 
ment of the charged nucleons inside the nucleus (convection current) and by the sum 
of the magnetic currents associated with the nucleons spin (magnetization current). 
The expressions of these operators are: 



P(r) = E^'5(r-r,), (22) 



fc=i 



for the charge operator. 



f(^) = E ^ ^ [^(r - r.)V. + V.5(r - r,)] , (23) 



k=i 2^fc 2 
for the convection current operator, and 



j"(r) - t i (/^^ + ^"^^) V X Kr - r.) .. (24) 

for the magnetization current operator. In the previous equations Mk indicates the 
rest mass of A;-th nucleon, /i^ and jj^ are the anomalous magnetic moment of the 
proton and the neutron respectively, ak is the Pauli spin matrix of the fc-th nucleon 
and r| = 1 or -1 according the k-ih nucleon to be a proton or a neutron. 

The expressions (0), (^3]) and (|2^ ) of the one-body nuclear charge and current 
are obtained as Fourier transform of a specific non relativistic reduction of a nucleon 
current of the type: 

J^\PfSf, P,S,) = u{Pf, Sf) j^ u{P„ S,) , (25) 



where the u are the nucleon Dirac spinors and P and S are respectively the four 
momentum and the spin of the nucleon. The procedure leading to the expressions 
(p2D-(p^) consists in retaining the zero and first order terms of an expansion in 
powers of q" = P^ - P^ . 

2.2 Coulomb distortion 

The possibility of making an experimental separation of the longitudinal and trans- 
verse response functions is very appealing because it allows one to investigate sep- 
arately the electron nucleus interaction and the nuclear excitation mechanism. The 
same nuclear model should be used to describe both responses, the only difference 
between them is due to the electromagnetic operator. 

We have seen that the Rosenbluth separation is possible if the expression ([T3|) 
of the cross section is valid. In this expression the kinematics is well separated from 
the nuclear structure since the PWBA has been used, in other words, since it has 
been assumed that the electron wave functions are plane wave solutions of the Dirac 
equation. In this approach the effects of the electromagnetic field of the nucleus 
on the motion of the electron are neglected. The validity of this approximation is 
quite good for few body systems or light nuclei, but it starts to become questionable 
already for nuclei as heavy as ^^C ||deF66| , |Cio80| , |Hei83|] . 



It is possible to solve the Dirac equation for an electron moving in a static 
Coulomb potential [Yen54|. This provides the exact solution for the elastic scatter- 



ing cross section. The situation is more difficult to handle in the case of inelastic 
scattering since the electromagnetic field generated by the nucleus is not any more 
static because it is produced by a transition from the ground to the excited state. 
In the most sophisticated treatment developed up to now, the one-photon exchange 
scattering amplitude (|^) is calculated using the electron wave functions obtained 
by solving the Dirac equation with a Coulomb potential generated by the target 
nucleus ground state. This treatment is called Distorted Wave Born Approxima- 
tion (DWBA) and its validity in the description of low lying states has been widely 
studied [[ilei83i . 

There are two major effects that the DWBA adds to the PWBA description. The 
first is produced by the attraction between the positively charged nucleus and the 
negatively charged electron. The asymptotically measured initial and final momenta 
ki and kf are smaller than the effective momenta of the electron around the nucleus. 
This is an interpretation of the phenomenon based upon a PWBA description of the 
process, since DWBA wave functions are not eigenstates of the momentum operator 
of the electron. In this simplified picture the incoming electron is accelerated by 
the nucleus while the scattered electron is decelerated. This produces a transferred 
momentum larger than that of the PWBA. An empirical expression of the effective 



momentum transfer is obtained in the eikonal approximation [Hei83 



/ Zahc\ , , 

geff = g 1 + 1.5^ (26) 
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Figure 2: Rosenbluth plot for ^^C(e,e')-'^^C and ^°Ca(e,e')^'^Ca reactions at the fixed 
excitation energy u = 300 MeV. We have indicated with um the Mott cross section. Note 
that the DWBA points are not perfectly aligned. 
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with i?eq = 1.24: A^'^. The difference between the PWBA momentum transfer and 
the effective one increases with the charge Z of the target nucleus and diminishes 
when the electron energy Ei increases. 

The second DWBA effect is of quantum mechanical nature and consists in the 
focusing of the electron wave function onto the nucleus causing an increase of the 
flux and an enhancement of the cross section. This effect is partially compensated 
by the reduction of the cross section due to the qes shift. 

In the DWBA expression of the cross section it is not possible to separate the 
nuclear from the electron variables. Furthermore there are interference terms be- 
tween longitudinal and transverse amplitudes. In a Rosenbluth plot the DWBA 
cross sections do not lie on a straight line. If the distortion effects are small, how- 
ever, it is possible to correct the experimental cross sections in order to have PWBA 
equivalent cross sections and to use them to make the Rosenbluth separation. 

Distortion effects on the QE excitation have been studied for ^^C and ^°Ca nuclei 
in ref. ||Co'87a||, where a full DWBA calculation has been done, and for ^°Ca and 



^°^Pb nuclei in ref. ||Tra88| , [Tra93|| , where an analytical expansion of the amplitudes 
up to (Za)'^ has been used. 

A first finding of these studies is that the DWBA cross sections without cor- 
rections are rather well aligned on a Rosenbluth plot. Within the experimental 
uncertainty their alignment is compatible with an exact straight line (see fig. ^. 
For this reason the alignment of the cross sections on a Rosenbluth plot cannot be 
considered a proof of the validity of PWBA. 

The major effect of the distortion on a Rosenbluth plot consists in a rotation of 
the line around the intersection point with the y axis. This does not modify very 
much the value of the longitudinal response but it has a rather big effect on the 
value of the transverse response function. These considerations are valid for a fixed 
value of the excitation energy. 



In ref. ||Tra88|| the full energy dependence of the responses has been studied. The 



distortion effects on ^°Ca consists mainly of a shift to higher energies and it can be 
rather well described in terms of gcff- The situation in 208p]^ jg j^iore complicated 
since the focusing of the wave function plays an important role. 

The effects of the distortion are small enough to be considered a correction to 
the PWBA, therefore it is reasonable to extract the responses from the cross section 
using the procedure outlined above. 



2.3 Nucleon electromagnetic form factors 

In the standard derivation, the nuclear current has been obtained under the as- 
sumption of pointlike nucleons. In the QE excitation the values of the energies and 
momenta involved are such that the nucleon internal structure plays a noticeable 

11 
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Figure 3: ^^C response functions. The bands show the uncertainty related to the choice 
of the the nucleon electromagnetic form factors. 
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role. This fact is taken into account by inserting nucleon electromagnetic form fac- 
tors in the expression of the current. For a nucleon of type r (either proton or 
neutron) with mass M , eq. (pSJ) becomes: 



J^^PfSf.PS^) =uiPf,Sf) 



77T/ 2\ M I • T-^2 W J fiu 



2M 



where 



a 



flU 



i^iYY-YY) 



uiP,,S^) (27) 



(28) 



and K^ is set equal to the anomalous part of the magnetic moment in units of nucleon 
magnetons (1.79 for proton and -1.91 for neutron). The terms F[ and FJ are the 
Dirac form factors and they are normalized such as: 



Ff(0) 

Fm 



1 





F|(0) 



1 
1. 



(29) 
(30) 



For a more direct comparison with the experiment it is convenient to use the Sachs 
form factors related to the Dirac ones by the expressions: 



Gg 



F[ 



G 



M 



k^f:^ . 



(31) 

(32) 



The effect of the nucleon finite size is implemented by means of a convolution 
of the pointlike currents (P^-(P^ with these nucleon form factors. Then, it is 
more convenient to work in q space where the nucleon form factors are simply 
multiplicative factors. 

The problem is now to find a functional dependence of the form factor in terms 
of q and u which is able to reproduce the elastic electron-nucleon data. Unfor- 
tunately the data are not selective enough to define a unique parameterization of 
this functional dependence. Different type of form factors can reproduce the elastic 
scattering data with the same degree of accuracy. 

The sensitivity of the nuclear responses to the different choices of the electro- 
magnetic nucleon form factors has been studied in ref. [|Ama93b|| . In that article we 
performed calculations of the QE responses with the same nuclear structure input 
but considering the nucleon form factors of refs. | Jan66 , Ber72| , [lac73| , [Hoe76 , SimSO 



We found that the momentum and energy domain of the QE peak lies in a region 
of strong variations of the nucleon electromagnetic form factors, and the changes of 
the nuclear responses can be as big as 12% at the peak values (see fig. ^. 

Some of the parametrizations we have used are old and have been ruled out 
by more modern measurements. Probably the 12% uncertainty we have quoted 
slightly overestimates the actual uncertainty on the form factor. In any case the 
study quoted above points out the need of specifying the nucleon form factor used 
when one compares the results of QE calculations with the experimental data. The 
ambiguities in the choice of the nucleon form factors set an intrinsic uncertainty 
threshold of the theoretical calculations. 
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2.4 Relativistic corrections 

We have already said that the expressions of the one-body operators (^)-([2^) are 



obtained by expanding the relativistic operator (|25| ) in powers of g, u and P, and 
retaining the zero and first order terms. 



In ref. |[Voy62|| it has been pointed out that a consistent treatment of both the 
charge and the current operators, up to the same order in the expansion, produces 
a correction factor for the charge operator. This factor is called the Darwin-Foldy 
term 

/o.^(l.^) (33) 

and it multiplies the electric Sachs form factor in the calculation of R^. This cor- 
rection reduces the longitudinal response as is shown in fig. |[ 

A new non relativistic expansion of the on-shell electromagnetic one-body current 
has been proposed in ref. |[Ama96a| for use in the region of the QE peak at high q 
values. The idea of this approach is to obtain expressions of the current operators 
which are not truncated in powers of q and u and hence applicable at high values 
of energy and momentum transfer. 

In a FG model the value of the momentum of a bound nucleon is at most of 
the order of the Fermi momentum P r^ k-p. Therefore for a typical value of k-p = 
250 MeV/c we have P/M < 1/4 and hence an expansion of the current in powers 
of P/M is justified. 

The details of the calculation can be found in refs. | Ama96a , Ama96b , Jes98|| 



and the final expressions of the charge and transverse current operators up to first 
order in P/M are: 

The first piece of the charge operator is the usual non relativistic contribution Ge 
multiplied by a factor larger than one. This factor produces (see fig. ^ an increase 
of the non relativistic result which goes in the opposite direction with respect to 
the traditional Darwin-Foldy correction, and it can be directly implemented as a 
multiplicative factor in the non-relativistic models of the reaction. The second piece 
is the charge of the spin-orbit operator, which is slightly different from the one 
arising in the usual expansions. Finally, the transverse current operator contains 
the two traditional pieces, the magnetization and convection current operators, but 
both are corrected by a multiplicative factor smaller than one. 

The results presented in fig. ^ show big relativistic effects at momentum transfer 
larger than 500 MeV/c. In these situations the value of the momentum of the emitted 
nucleon is comparable with its rest mass M , therefore a relativistic treatment of the 
motion of the emitted nucleon is necessary. This involves the nuclear structure part 
of the problem and we shall discuss this point more in detail in section p.2| . 
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model. The dashed lines present the result obtained by using the usual non relativistic one- 
body operators. The dotted curves include the Darwin-Foldy term (only in the longitudinal 
response). The full lines have been obtained considering the relativistic corrections to the 
operators. 
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2.5 Two-body currents 

We have treated the problem of the electron excitation considering the nuclear cur- 
rent produced only by individual nucleons. This limitation is inconsistent already at 
the level of the standard derivation. Inserting the expression (^ into the continuity 
equation (^ one obtains: 

V-J(r) = -i[H,p\ = -^\T + V,p] , (36) 

where T is the kinetic energy and V the interaction term. 



It can be shown that by using the expressions ( P^ and ( P5D for the charge and 
convection current operators the equation 

V-f(r) = -t[T,p] (37) 

is satisfied. It is obvious that a new term must be added to the one-body current 
in order to verify the continuity equation with the full hamiltonian. This new term, 
J"""^, must satisfy the equation 

V-J™(r) = -z[V,p] (38) 

and it is a two-body operator. 

The continuity equation shows the theoretical needs for including two-body cur- 
rents, but it does not define them in an unique way, because all the operators which 



are divergenceless, for example the one-body magnetization current (|2^) , are not 
restricted by this equation. 

Since there is not a unique way to derive these two-body current operators, dif- 
ferent approaches have been proposed in the literature. In any case, it is commonly 
accepted that the most important terms are those related to the exchange of mesons 
and more specifically to the exchange of a single pion. 

We have calculated the contribution of the Meson Exchange Currents (MEC) 
following the approach of refs. ||Gar76|, [b'ri77|, |Van81||. Using the Feynman rules we 



have evaluated the contributions of the three diagrams of fig. ^ which we call seagull 
(a), pionic (b) and A-isobar (c) terms. 

For the evaluation of the MEC terms we need a model for the pion-nucleon, 
photon-nucleon- A and pion-nucleon-A interactions. The models of the pion-nucleon 
interaction commonly used (pseudo-scalar and pseudo- vector couplings) produce the 
same result in the non-relativistic limit. 

We define a function h{r) as the Fourier transform of the dynamical pion prop- 
agator 

u \ f ^'^ i^.N(fc,e)e-'-('-') 
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Figure 5: Feynman diagrams of the three MFC operators we consider: (a) seagull, (b) 
pionic and (c) A-isobar. 



where we have indicated with F^n the pion-nucleon form factor and with e = {A.E)i 
the energy of the exchanged pion obtained as the difference between the energies of 
the final and initial single states of the l-th nucleon. 

In momentum space we express the seagull and pionic currents respectively as: 
f{q, u) = 47r^ Fs{q, u) J] [r^ x n], e'^^'^ak en ■ VkK^u - n) (40) 



and 



/. 



r(q,a;) = -A7C^F^{q,uj) Y. i^k x ri]^ U\ e'^'^ 

TTT'^ I, ; — 1 J 



k,l = l 



(Tk ■ V/i(r - rfc) V [ai ■ V/i(r - n)] 



(41) 



The situation for the A-isobar current is not so well defined because formulations 
based upon static quark models ||Che71| , [Hoc73| , |Loc75| , |Ris79|| or chiral lagrangian 
[Fec69| , |Van81| , |Alb90| | give different expressions. We have adopted the point of view 
of the first group of authors and use the following expression: 



i^{(i,r) 



-iC^F^{q,u) J2 e'^^^'q 



k,l = l 

k^l 



X |[rfe X Ti]3 (Tk X VkO-i ■ Vkh{vk - r^) - ArfVkCri ■ ^kh{vk - r;)| . 

In the above equations /,r = 0.079 is the effective pion-nucleon coupling constant, 
m^ is the pion mass and 



CA = 47r 



J n 



ml 25M(Ma - M) 



(42) 



with Ma the A mass. Finally, Fg, F^^, Fa show the dependence on the electromag- 
netic nucleon, delta and pion form factors. As we have already discussed there are 
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q [MeV/c] cu [MeV] AF, 



s 



300 
400 
500 



50 
90 
100 



15.7 
25.6 
40.2 



AF. 



5.3 
8.3 
12.3 



Table 1: Percentile uncertainty on the value of the peak of the response produced by the 
different choices of form factors on seagull and pionic currents |Ama93a]. 



theoretical ambiguities in the choice of these form factors. In table |I| we show the 
effect of these ambiguities on the values of the peak of the MEC responses ||Ama93a[l . 



For the seagull current, the possible functional dependence of the F^ form factor, 
generate an appreciable uncertainty at the peak energy. For q = 550 MeV/c, the 
uncertainty on the maximum is ~ 40% of its value. A similar situation is found 
for the electromagnetic pion form factor F^ included in the pionic current even if 
the uncertainty is considerably smaller than in the seagull case. Table shows for 
example that, at g = 550 MeV/c, the uncertainty of the peak value is ~ 12%. 

For consistency with the one-body currents, we use in our calculations the fol- 
lowing expressions for F^ and F^^: 



F^ 



G^, - G 



E5 



7r7 



1 + (g^ — iu'^)/Tn'l 



(43) 
(44) 



where rup is the mass of the p-meson. 



The situation is even more complicated for the A current since the electro- 
magnetic form factor F^, and the constant C^ are model dependent. The major 
uncertainty is related to Ca, but a discussion of this problem is beyond the aim of 
the present work. The expression (^21) of Ca we have chosen is widely used in the 
literature [ Hoc73 , K.is79 , K.is84 , [Sch89|] . For the form factor Fa, following the static 
quark model, we use: 

Fa = 2G^ = Gyi — G]y[ . (45) 



Finally, we would like to comment on the pion-nucleon form factor 



F^n{k,€] 



A' 



TTIZ 



A2 + A;2 - £2 



(46) 



present in the pion propagator of eq. (^). We have verified [|Ama93ai |Ama93b( 
that, in the QE peak region, for the values of A commonly accepted (~ 1 GeV), the 
results are very close to those obtained considering simply F^n = 1, which is the 
value we have adopted. 

The MEC we have presented above act only on the transverse response. In 
principle MEC contributions are present also in the longitudinal response, but it 
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Figure 6: Ip-lh MEC contributions to the excitation of the Ca transverse response at 
g=550 MeV/c. The full line shows the interference term between one-body current and 
MEC. The dashed line represents the pure MEC contribution. 



has been shown |P^al97|| that their effects are very small, of the order of 3% at the 
energy of the peak. For this reason in the following we shall neglect the corrections 
to the longitudinal response and we shall discuss only the effect of MEC on the 
transverse one. 

The MEC are two-body operators, therefore they can produce the excitation of 
Ip-lh or 2p-2h pairs. In the excitation of Ip-lh pairs, the transition amplitudes of 
MEC and those of one-body currents interfere. These interference terms produce 
the larger contribution of the MEC to the Ip-lh response. 

In fig. ^ we show the contribution of the MEC to the transverse response of 
^°Ca for g=550 MeV/c. While the contributions of pure MEC are positive for 
all excitation energies, the contribution of the interference term between one-body 
current and MEC is positive at small energies and become negative at energies above 
80 MeV. 



In agreement with ||Koh81|| we find that the effect of the pionic and seagull 
currents on the responses have opposite sign. The pionic current diminishes the 
pure one-body response, while the seagull increases it. In fig. |^ we show this fact for 
the ''^Ca nucleus, but we found analogous results also for the ^^C response ||Ama93a 



Considering only seagull and pionic currents we obtain an increase of the transverse 
response in the peak region ||Ama92|, |Ama93b|l. The inclusion of the A isobar current 
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body responses. The remaining Hnes show the results obtained by adding the vari- 
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Figure 8: 2p-2h responses of ^^C for q = 550 MeV/c. Panel (a): Seagull plus pionic 
contributions. The dashed line is the term corresponding to the two particle emission 
while the dotted line has been obtained for the case of one particle in the continuum 
and the other one on an excited bound state. The full curve represents the total 2p- 
2h response. Panel (b): Two nucleon emission contributions induced by the A isobar 
current. The dashed-dotted line corresponds to the emission of a proton-neutron pair. 
The dotted and dashed lines give the contribution of the two-proton and two-neutron 
emission, respectively, and are overlapping in the scale of the figure. The solid line is the 
total response. 

gives a negative contribution. The final results consists in a lowering of the one- 
body responses. The effect of the seagull and pionic currents becomes smaller as 
q increases while that of the A isobar current remains roughly constant. For this 
reason, the total lowering effect of the MEG increases its magnitude with q. 

In all the nuclear models we have used to calculate the response, the nuclear 
Ip-lh final states are orthogonal to the 2p-2h final states. This means that in our 
approach there is not interference between Ip-lh and 2p-2h excitation amplitudes 
and, as a consequence, the 2p-2h responses add to the Ip-lh ones. 

The two body MEG operators can lead to different kinds of 2p-2h final states. 
A first classification concerns the number of particles which are in the continuum. 
A first type of 2p-2h state is produced when none of the two particles is in the 
continuum. These are low energy excitation and they do not contribute in the 
QE region. A second kind of excited state is formed when one particle is in the 
continuum and the other one is lying on a bound single particle state. Finally the 
last type of excited 2p-2h state is the one which has both particles in the continuum. 

In the panel (a) of fig. § we show the behavior of the 2p-2h responses for the 
various types of final states. In this calculation only the seagull and pionic currents 
have been used. The contribution of the response where only one particle is in the 
continuum rapidly reaches a plateau value and it remains roughly constant with the 
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total responses including Ip-lh and 2p-2h final states. The dotted lines give the one-body 
response and the dashed lines the 2p-2h contributions. The dashed-dotted lines have been 
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increase of the excitation energy. Since the number of particle-hole configurations 
entering in the excitation of this part of the response remains constant even if the 
energy increases. The two-nucleon emission contribution behaves quite differently: 
its value increases steadily with the value of the excitation energy. This is related 
to the enlarging of the phase space available for the final states of this reaction. 

The other possible classification of 2p-2h final states is related to the kind of 
particles which are emitted. Seagull and pionic currents are generated by the ex- 
change of charged pions. The A isobar current is active also when a neutral pion 
is exchanged. Then, while the seagull and pionic currents can lead to final 2p-2h 
states where only a proton and a neutron pair is emitted, the A isobar current can 
emit also a pair of like nucleons. 

The contribution of these different kind of excited states on the 2p-2h response 
of the A isobar current is shown in panel (b) of fig. ^. The two-proton and two- 
neutron emission responses are indistinguishable at the scale of the figure. This is 
due to the fact that ^^C nucleus is practically symmetric in terms of proton- neutron 
single-particle structure. The number of protons and neutrons is the same, and the 
small differences in the single-particle energies and in the wave functions affect very 
little the response. The dominant emission channel remains the proton-neutron one, 
and its relative contribution becomes more important with increasing energies. 

In fig. 1^ we present the total QE transverse responses for ^^C and ^°Ca. The 
2p-2h responses are rather small below the peak but become important at higher 
energies. 

A detailed analysis of the MEC contribution |[Ama94a| shows that the pionic 



and seagull currents produce an increase of the values of the responses in the peak. 
The inclusion of the A isobar current has an opposite effect. If for g=300 MeV/c 
the total responses are similar to the one-body responses, for higher values of q they 
are even smaller. There is not big difference between the two nuclei, though in ^'^Ca 
the modifications are slightly bigger. The greater effect we found is less than 5%, 
certainly within the range of uncertainty related to the choice of the nucleon and 
mesons form factors. 



3 The nuclear excitation 

In this section we shall present some of the models commonly used to describe the 
nuclear QE excitation. In this region the excitation energy of the nucleus is well 
above the continuum threshold, therefore at least one particle is emitted from the 
nucleus. Due to the relatively high values of energy and momentum transfer the 
QE excitation is not of collective type like vibrations or rotations, but the process is 
dominated by the single particle dynamics. The MP model is the obvious starting 
point of any description of this excitation. We shall first present the basic ideas 
of the MP model and in subsequent subsections we shall discuss the validity of the 
approximations made in this treatment. 
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3.1 Mean field models 

In the MF model the many-body nuclear states are Slater determinants of sp wave 
functions. In the ground state all the sp states below the Fermi level are occupied 
and all the states above it are empty. The excited states of the system are generated 
by promoting nucleons from below to above the Fermi level. These basic ideas of 
the MF approach can be used with different set of sp bases, and this differentiate 
the various MF models. In our presentation we shall treat the shell model where 
the sp basis is generated considering the finite size of the nucleus and the FG model 
where the system is considered to be infinite and translationally invariant. 

3.1.1 The Continuum Shell Model 

The basis commonly used to describe finite nuclear systems is composed of sp states 
characterized by the orbital angular momentum /, the total angular momentum j, 
its projection m on the z axis and the isospin z axis projection, t. For discrete sp 
states it is also necessary to specify the principal quantum number n, 

\nljmt) = <,,v„(r) = Rii^{r) yijmij-) Xt (47) 

where 

yiU^) = E {if^UM M^) xs ■ (48) 

flS ^ 

Here {lfi^s\jm) is a Clebsh-Gordan coefficient, Xs is the spin function, with s the 
spin z axis projection, and Xt is the isospin function. For continuum sp states their 
energy e must be given 

\eljmt) = ip'ij^ir, e) = Rjj (r, e) yijm Xt ■ (49) 

One possibility of generating the set of sp states is to use a Hartree-Fock pro- 
cedure as is done in ref. ||Koh83|| . This implies the choice of an effective nucleon- 



nucleon interaction and the solution of a system of non linear Schrodinger-like dif- 
ferential equations. 

We use a different approach called Continuum Shell Model (CSM). We fix a 
central real potential V{r) acting on all the nucleons, and we generate the sp states 
( ^7D and (PPP by solving the traditional radial sp Schrodinger equation with the 
appropriate boundary conditions on the radial parts. The solutions of this equation 
should be regular at the origin. At high r values bound states have exponential 
decay, while continuum states have oscillating behavior. 

As we deal with one- and two-body operators, the only possible final states are 
of Ip-lh and 2p-2h types, with the particles in the continuum, which in this model 
are, respectively, given by: 

\f3,JMj) = \{eplpjp),{nhlh3hy,JMj), (50) 
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Figure 10: Sensitivity of the responses to the choice of the single particle wave func- 
tions. Full lines represent the results obtained with the Woods-Saxon set of refs. 
| Co'84 , po'85 , po'87b|| . Dashed lines correspond to a new parameterization aiming to 
reproduce the charge density distributions. Thick full lines represents the empirical 
charge density distributions reconstructed using the Fourier-Bessel parameters given in 
ref. fpeJS?!] . 
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\P, JMj) = I [(epjpjp,), (epjpjp^); Jp] , [(nhjh^jhi), (rihjhijhi); A] ; JMj) . (51) 



In the case of one-body operators, only the Ip-lh final states can be excited. The 
sum over (/?J) in eqs. (|20|) and (|2l|) goes on all the possible Ip-lh configurations 
allowed by the angular niomentuni and parity composition rules. This means that 
\jp—jh\ < </ < jp+jh and (— l)^p+^'' should be equal to the parity of the corresponding 
electromagnetic transition. The transition amplitudes can be written as: 



OB( 



(/?, J II Oy^(g) II 0) = {tplpjp II Oy^(g) II n^hjh) S{uj - e^ + e,) . 



(52) 



In case of two-body operators the expressions are slightly more complicated, 
and it is possible to excite both Ip-lh and 2p-2h final states. In the first case the 
same angular momentum, parity and energy rules as before apply and the transition 
amplitudes can be written as: 



(/3, J II 0™(g) II 0) = E E (-1)^^+^^+^+^ V(2Ji + l)(2J2 + l) I 

a<ep J1J2 { 



Jl J2 J 

ep J,J2 K ^^ ^P -^^ 

{{nJo,ja),{eplpjp)]Ji II 0™(g) || {uhlhjh), {riaLja); J2) + 
{-iy-+^'^+J^+^{{nJ^j^),{eplpjpy,Ji II 0™(g) || (n^J.), (n^4jh); ^2) ^ 
6{uj -ep + eh) , (53) 

where we have used a 6j coefficient. 

If the final state has two particles in the continuum the transition amplitude is: 
{(3, J II 0™(g) II 0) = (-l)M+...+^P S{uj - ep, - e^, + e,, + e,J 

((ep2^P2Jp2)> (epi^PiJpi); Jp II 0]^{q) \\ (nhjhjh^), {nh^h.Jh^); Jh) + (54) 

^_iyh^+3h^+Jh+i {{epjpjp,,), {epjpjp,); Jp \\ Oj^(g) || {ntjh^jhi), {ntjhjh^); Jh) 



where now \J„ 



Jh\ ^ J ^ Jp + Jh- 



The above expressions show that in the MF model, the complicated problem 
of evaluating many-body transition matrix elements has been reduced to the much 
simple problem of calculating sums of one- and two-body transition matrix elements. 
The calculation of the sp transition matrix elements is long and tedious, but it does 
not present fundamental difficulties and it can be carried on without any further 
hypothesis or approximation. The explicit expressions of the matrix elements can 
be found in ref. ||Ama93a[ . 



We have studied the sensitivity of the CSM results to the choice of the potential 
parameters |[Ama934 |Ama93 b| ] , the only input of the model. We have performed 
CSM calculations with two sets of parameters. The first one taken from the literature 
has been fixed to reproduce the sp energies around the Fermi level and it has been 
used to describe low-lying and giant resonance excitations ||Co'84|, |Co'85|, |Co'87b| . 
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The second set of parameters has been fixed to have good description of the enipir- 

The comparison of the results 
The difference between the 



ical charge density distribution |[Ama93a| , |Ama93b(] 
obtained for one-body responses is shown in fig. 



10 



results obtained with the two parametrizations becomes smaller at high values of q 
and u. This indicates that in the QE peak the surface effects are not important. 
The differences between the two results are compatible with the uncertainty related 
to the choice of the nucleon form factors. All the CSM results we shall present have 



been obtained with the set of parameters of refs. 



rVHTE 



3.1.2 Fermi gas model 



An alternative MF model often used ||Czy63| , |Mon69| , [Koh81| , [Van81| , [Bou91| , |Gar92| , 
Gil97|| to describe the QE excitation is the FG model. In this model the nucleus 
is treated as an infinite system of non-interacting fermions. Both the volume Q of 
the system and number of particles A are infinite, but the density p = Q/A remains 
finite. In this model the meaningful quantities are the responses per nucleon. 

Since the nuclear system is translationally invariant, the adequate set of sp wave 
functions is formed by eigenstates of the momentum 



\kst) 




(55) 



In the nuclear ground state all the sp states with k < kp are occupied. The Fermi 
momentum kp is related to the constant density of the system by the expression: 



3 2 ^ '/' 

-TT p 
2 ^ 



(56) 



from which, assuming a hard sphere model of the nucleus, p = (|7rrQ j and for tq 
= 1.12 fm, one obtains the well known values kp = 1.36 fm~^ or 268 MeV/c. 

In analogy with the CSM the excited states are generated by promoting a particle 
from a level below to above the Fermi level. While in the CSM model the angular 
momentum of the system should be conserved, in the FG model the total momentum 
of the system is conserved. This means that the relation q + Pi = P/, where p, and 
Pf are respectively the initial and final nucleon momenta, should always be satisfied. 
For the Pauli exclusion principle |p/| should be greater than kp, and therefore those 
excitations such that |q + Pj| < kp are forbidden. If |q| > 2kp there are no forbidden 
excitations, i.e. every nucleon can be emitted to the continuum. 

In the FG model there are not discrete sp states, therefore it cannot be used 
to study low-lying nuclear excited states. Furthermore the excitation spectrum of 
this model does not predict the giant resonances, since they are collective surface 
vibrations of the nucleus in the continuum and the FG does not have surface. The 
QE region is however well suited to be studied with the FG model. The excitation 
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energy is such that the nuclear finite state is always in the continuum and the 
nuclear dynamics is of single-particle type which one expects to be more related to 
the volume than surface effects. 

The expression (^) of the sp states allows one to perform a large part of the cal- 
culation of the responses in an analytic way. The final expressions for the responses 
to be handled numerically are much simpler than those obtained in the CSM. These 
expressions can be found in 



Ama93a, A.ma94b 



The shape of Rl and Rt changes when q become larger that 2kp. For q < 2kp 
the response as a function of the excitation energy u has first a linear increase then 
it decreases as a parabola. For q > 2kp the response has a parabolic behavior with 
a maximum in tu = g^/2M. The energy and momentum conservation implies that 
for a fixed value of q there is a maximum value of u where the response is different 
from zero. Of course the experimental data do not show these behaviors. 

We have studied the reliability of this model in the description of the QE excita- 
tion by comparing the nuclear responses obtained with CSM calculations with those 
obtained with the FG model. The only input in the FG calculations is the value of 
kp, therefore of the nuclear density. The results we have obtained for the traditional 
nuclear matter value of kp, are represented by the dotted lines of fig. 0. These 
responses are quite different from those calculated within the CSM (full lines). The 
maxima are smaller and the widths are broader. It is remarkable that the positions 
of the peaks are roughly the same for both calculations. 

The value of kp used in these calculations has been obtained under the assump- 
tion that the density of the FG system corresponds to the central density of the 
nucleus. A better approximation consists in considering an average nuclear density 
by using a value of the Fermi momentum evaluated as: 

(kF) = -vrM L (57) 

In our calculations the effective momenta for ^^C and ^'^Ca have been obtain using 
the proton density produced by the Woods-Saxon potential (see fig. |T0|). We found 
(A;f)=215 MeV/c for ^^c and (A:f)=235 MeV/c for ^^Ca. The use of the effective 
Fermi momentum greatly improves the agreement with the CSM results as is shown 



in fig. |TT|. As expected the agreement is better the heavier the nucleus is and 
the higher the momentum transfer becomes. This is a further indication that, under 
these conditions, the QE excitation is mainly a volume effect and the nuclear surface 
plays a minor role. 



The Local Density Approximation (LDA) is a prescription commonly used | |Gil97| 



to correct the FG calculations for the fact that the nuclear density is not constant. 
The philosophy underlying the LDA is that, at any point of the finite nuclear volume. 
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Figure 11: Comparison between shell model responses in ^^C (full lines) and FG responses. 
The dotted lines have been calculated using the Fermi momentum value extracted from the 
nuclear central density A;f=268 MeV/c, the dashed lines with A;f=215 MeV/c calculated 
with eq. (|57[) and the dashed-dotted lines with the LDA, eq. (p3). 
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Figure 12: Interference terms between pionic and seagull MEC and one-body currents 
for 12C (left panels) and 40Ca (right panels). The full lines have been calculated within 
the shell model, while the dashed lines show the FG results obtained using the average 
values /cf=215 MeV/c and A;f=235 MeV/c for the two nuclei, respectively. The pionic 
contribution is negative and that of the seagull current positive. 
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the response can be described by the nuclear matter response calculated at the 
corresponding density. Using eq. (p6]) we evaluated the LDA responses as: 



^"°^(g,^) = / 



3 i?^G(g,a;,fcp(r)) 



dV 
2 



SttM 



d\k\r)R^^{q,u},kF{r)) 



(5^ 



where we have indicated with R^'^ the FG response, and with Q{r) the nuclear 
volume defined as Q{r) = A/p{r). 

The dashed-dotted lines of fig. |Tl] are the LDA responses calculated using the 
shell model proton density. The shapes of the LDA responses are more similar to 
those of the GSM than those obtained with the FG. While these latter have a sharp 
cut in energy, the LDA produces a tail, as in the GSM calculations. In spite of the 
qualitative similarity between LDA and GSM responses, from the quantitative point 
of view the agreement between these two results is rather poor. 

We have discussed so far the FG results obtained for one-body operators only. 
We have also calculated the MEG contribution in this model, by considering the 
most important terms, i.e. those produced by the interference with the one-body 
currents. In fig. |12| we compare the interference terms of pionic and seagull currents 
with the one-body currents calculated in the GSM and in the FG model (using 
effective momenta). Also in this case the FG gas model produces results very similar 
to those obtained in GSM, especially for high values of q. 



3.2 Relativist ic corrections 



In section p.4| we discussed the relativistic corrections to the expressions of the one- 
body electromagnetic operators. The results presented there have been obtained 
using the relativistic energy- momentum relation E'^ — p'^ = M^. The GSM and the 
FG model are non relativistic, therefore the use of the expressions (0) and (|35| ) is 
not correct within these models. 

A prescription to consistently implement the expressions ( plD and (|35|) within a 
non relativistic description of the nucleus has been presented in ref. | |Ama96a |. The 
basic idea consists in calculating the wave function of the emitted particle of energy 
Cp = e/i + (^ by inserting in the Schrodinger equation the relativistic version of the 
kinetic energy operator. This is equivalent to solving the Klein-Gordon equation 
instead of the Schrodinger one. The asymptotic momentum of the ejected nucleon 
is given by: 



P 



(ep + Mf -M^ = ep{ep + 2M) = 2M 



uo 



' + 2m) + ^4' + -M 



eft 



(59) 



In the non-relativistic regime uj is much smaller than the nucleon mass, therefore 
the 1/M terms can be neglected and the non relativistic relation p^ = 2M{uj + th) 
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Figure 13: The full lines show the results of the exact relativistic FG model. The dotted 
lines have been obtained considering the relativistic corrections to the one-body operators 
together with the relativistic kinematics describing the motion of the emitted nucleon. 
The dashed lines show the non relativistic FG results. 
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is recovered. A better approximation is obtained by neglecting only the t^ 

term. This approximation can be implemented in a non relativistic model by making 
the substitution lo — >■ a;(l + a;/2M). This modification of lo should be used only for 
the calculations of the the responses |[Alb9CI||, but not in the nucleon form factors, 



because the energy and momentum transfer are still q and uj (the only difference 
is that the momentum of the ejected particle is calculated by using the relativistic 
energy- momentum relation). 

It is possible to test the validity of this prescription by comparing the results 
obtained with those of the relativistic FG model. In this last model the sp wave 
functions are free solutions of the Dirac equation and the relativistic expressions of 
the current operators can be used without doing any approximation in the calcula- 
tions of the responses. We show in fig. |T^ this comparison and we see that, even 
at high values of q where the relativistic effects are more important, the agreement 
between the results obtained with the prescription presented above and those of 
the relativistic FG is excellent. The prescription can also be used in CMS calcula- 
tions and, for high q where surface effects are negligible, good agreement with the 
relativistic FG model results has been found ||Ama96a| . 



This procedure for including relativistic effects has been successfully applied 
to the study (e,e'p) processes from polarized deuterium ||Jes98|| and heavier nuclei 



A-ma96b| , |Ama98a[] . More recently the same expansion procedure has been applied 



to the MEG, although in this case the complexity of these current operators make 
the calculations more involved. The first results ||Ama98c| indicate a simple way 



of correcting the MEG for relativistic effects by using meanly multiplicative factors 
which depend on q and lo. 

From the above investigation it turns out that the relativistic kinematics produce 
a reduction of the width of the inclusive one-body responses and a shift of the 
peak positions towards smaller values of tu. In fact, the peak position is located 
at cu = (g^ — a'^)/2M, which for cj -C M can be approximated by a; = g^/2M, a 
larger value. The reduction of the width of the QE responses is noticeable in the 
upper limit of the lo region already for q = 500 MeV. For instance, for a value of 
kp = 235 MeV, about the average density of ^'^Ga the maximum allowed value of uj 
is ~ 225 MeV in the relativistic FG, against ~ 260 MeV in the non relativistic FG 



see ref. |[Ama96a| for details). 



In summary, the relativistic corrections to the electromagnetic nuclear responses 
are important for high q and the model developed for "relativizing" the currents can 
account for these corrections in a way which is easy to implement in already existing 
non-relativistic models. 



3.3 Random Phase Approximation 

In the MF model the many-body hamiltonian is a sum of single particle hamiltonians. 
This assumption neglects part of the interaction between nucleons which, for this 
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reason, is called residual interaction. Under certain approximations, the effects of 
the residual interaction on the description of the exited states of the many-body 
systems are considered within the Random Phase Approximation (RPA). 

The basic hypothesis of the RPA is that the nuclear excited state |/) can be 
described as a linear combination of particle-hole (p-h) and hole-particle (h-p) exci- 
tation of the ground state \i): 

\f) = Y1 {^IhO't^h - Y^y^alap) \i) , (60) 

ph 

where a^ and a are the creation and annihilation operators and p (h) indicates a 
particle (hole) state. The RPA amplitudes X and Y are provided by the theory and 
are normalized such as: 

E [ixi,r - iYj,r] = 1. (61) 

ph 

Once the X and Y RPA amplitudes are known, the expression ( pDj ) is used in 
the eqs. (|15D and (0) to calculate the transition amplitudes and the responses. In 
analogy to what happens in the shell model case, the evaluation of the complicated 
many-body transition matrix element is reduced to the calculation of a set of one- 
body transition matrix elements whose contributions for a specific transition are 
weighted by X and Y. 

The RPA has been implemented in nuclear physics in various ways, but for the 
application to the QE region it is necessary to consider excitations to the continuum. 
This means that in eq. ( |BUD the sum on p should be understood as a sum on the 
Ipjp and tp quantum numbers characterizing the particle state and an integral on its 
energy ep. 

The RPA can be applied also to an infinite system of fermions, and in fact it is 



in such a system that the RPA was first formulated ||Boh53|| . Whence the quantum 
numbers characterizing the sp states are the momenta. Thus the two sums of eq. 
(BUD become three dimensional integrals on the momenta. 



Following the approach of Section |3]2 we shall discuss separately the RPA results 
in finite and infinite systems. 

3.3.1 Finite systems 

The various implementations of the RPA to study nuclear excitations in the contin- 
uum, the Continuum RPA (CRPA), can be classified in two different approaches. 

In a first approach, the nucleon-nucleon effective interaction, used first in a 
Hartree-Fock calculation to generate the set of sp wave functions, is also used to 
perform the CRPA calculation. This approach, called self-consistent, is usually 
based upon the zero-range Skyrme interaction ||Sky56| ]. In QE | ]Cav84| , |Cav9(J| , |Van95[ | 



calculations various parametrizations of this interaction have been used. All of them 
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reproduce quite well the ground state properties of the doubly closed shell nuclei, 
and give a reasonable description of the giant resonances [|Vau72| , |War87| , ^ar93 . 

The second approach is based upon the Landau-Migdal theory of finite Fermi 
systems ||Mig57|| . This theory assumes a good knowledge of the ground state of the 
system and the RPA equations describe density fluctuations around the equilibrium. 
In the case of nuclei the ground state is usually described by means of a Woods- 
Saxon potential. Within this approach various kind of interactions have been used: 
from the traditional Landau-Migdal interaction of zero-range type up to finite range 



interactions [ Bri87 ] like the polarization potential [|Co'88 | and the G- matrix | Bub91 , 



Jes94| 



In fig. |14| we show a comparison between MF longitudinal responses (full lines) 
and CRPA responses calculated with the Landau-Migdal interaction (dotted lines). 
The peaks of the CRPA responses are about 10% smaller with respect to the peaks of 
the MF responses. Furthermore their position is moved towards higher values of the 
excitation energy. This effect is more pronounced at high q values. Similar effects 
have been found also in self-consistent calculations with the Skyrme interaction 
Dav8|, |Cm^ , |Van95|| . 



The applicability of zero-range interactions in the QE region is quite doubtful 
because of the high values of q and uj. The zero-range interactions are constant 
in momentum space, while finite-range interactions go to zero at high q. For this 
reason calculations done with zero-range interactions overestimate the RPA effects 
in excited states at high q and uo. The implementation of finite-range interactions 
in CRPA calculation is quite complicated because of the need of including exchange 
diagrams which cannot be easily handled in RPA. 

The polarization potential Pines et al. 



Pin88a, Pin88b] has the nice feature of 



being a finite-range interaction constructed to be used in RPA calculations with only 
direct terms. The effect of the exchange terms is considered in an effective man- 
ner with an appropriate choice of the parameters of the interaction. The dashed 
lines of fig. |14| show CRPA results obtained with the polarization potential ||Co'88| 



As expected the RPA effects are smaller than those obtained with the zero-range 
interaction and they become negligible with increasing values of q. More elabo- 
rated CRPA calculations which consider finite-range interactions in both direct and 
exchange terms |Pel85| , |Bri87| , ^hi89| , Pub91| , |Jes98|| produce analogous results. 

In spite of the differences in the technology of the solution of the CRPA equations 
and in the input, the various calculations in the QE region show that the RPA 
correlations lower the longitudinal responses and enhance the transverse ones (see 



right panels in fig. |l^). When a finite-range interaction is used, the RPA corrections 
to the MF results are quite small, within the uncertainty band related to the choice 
of the nucleon form factors. 



3.3.2 Infinite systems 



In infinite systems the QE excitations are often calculated ||Alb84| , |Alb87aj , |Alb87b| , 
Alb89| , |Alb9q , PeP93| , 1^11971 using the so-called Ring Approximation (RA). The RA 
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Figure 14: Responses calculated in CSM (full lines) and with Continumm RPA. The dot- 
ted lines have been obtained with a zero range Landau-Midgal interaction | Co'84 , Co'85f| 
while the dashed lines have been obtained using the finite range polarization potential 



| Pin88a| , piii88b|] . 
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is a RPA where the exchange diagrams are neglected. This allows one to obtain very 
compact expressions. In RA, the polarization propagator 11, whose imaginary part 
is proportional to the response, can be calculated as: 

where Vph is the effective nucleon-nucleon interaction and 11° the free propagator. 

Since the RA calculations in an infinite system are rather easy, it has been 

commonly accepted that RPA results can be simulated by RA by properly changing 

the parameters of the effective interaction. This issue has been investigated in refs. 

Bau98a, Bau98b| where RPA and RA transverse responses have been compared. 



The RPA responses have been calculated using a residual interaction containing 
a finite range term produced by the exchange of the vr and p meson plus a zero- 
range term of Landau-Migdal type. The RA calculations have been done using the 
same interaction but modifying the spin-isospin term of the zero range part, the Qq 
parameter, in order to best approximate the RPA results. The comparison of the 
two calculations has been done for various values of the momentum transfer q. The 
values of g^ needed to reproduce the peak positions of the RPA responses (squares) 
and those needed to reproduce the non energy weighted sum rule (triangles) are 
shown in fig. 0. Only around q =400 MeV/c it is possible to find a value of g^ 
able to reproduce both properties, but in general the results of the two theories are 
incompatible. 

The choice of the effective nucleon-nucleon interaction remains a major problem 
even in full RPA calculations for infinite systems. In finite system RPA calcula- 
tions the interaction is fixed to reproduce properties of low-lying excited states and, 
within the same RPA approach, it is afterwards used to calculate the QE response. 
This procedure cannot be used in infinite systems and therefore the interaction is 
taken from finite systems RPA calculations. From the theoretical point of view 
this procedure is inconsistent since the effective interaction is properly defined only 
within the framework of the effective theory where it is used. Usually the situa- 
tion is even worse since the interaction taken from finite systems RPA calculations 
is used to perform RA calculations in infinite systems. A detailed study of these 
inconsistencies has been done in refs. ||Bau98a], [Bau98b . 



In spite of these theoretical inconsistencies the RA results show effects rather 
similar to those of the CRPA calculations. The reason of this agreement is due to 
the small influence of the effect investigated especially at high q values. 



3.4 Final state interaction 

In any nuclear disintegration process the emitted nucleons interact with the residual 
nucleus. This obvious physical effect, which takes the name of Final State Interaction 
(FSI), is not considered in a FG model where the wave function of the emitted 
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Figure 15: Momentum dependence of the values of the parameter g^ to be used in RA cal- 
culations in order to reproduce the peak positions (squares) and the non-energy weighted 
sum rule (triangles) corresponding to the RPA responses. The dotted line gives the value 
g'o = 0.76 used in the RPA calculation. The dashed-dotted line shows the value q'q = 0.717 
which reproduces some low-energy properties by means of a RA calculation. 



nucleoli is a plane wave. In the CSM the emitted particles interact with the rest 
nucleus through the average nuclear potential. Strictly speaking this is part of the 
FSI. The CRPA calculations consider the mixing of the various nucleon emission 
channels and this is another part of the FSI. In the literature, the name FSI has 
been used to indicate those effects on the final state of the system which are beyond 
the MF and RPA descriptions. In this section we shall discuss them. 

A microscopic description of the FSI should be done with theories beyond the 
RPA. One of these theories is the second RPA | Yan82| , [Yan83 , Pro90|| which describes 
the nuclear excited state as a linear combination of Ip-lh and 2p-2h excitations. In 
ref. [pro87|] second RPA calculations in the QE region have been performed in 
the ^^C nucleus for low values of the momentum transfer. In these calculations 
the continuum has been discretized. The excessive numerical effort needed imposed 
a limitation of the configuration space used, which turned out to be too small to 
provide a realistic description of the QE excitation. In spite of this limitation, the 
results obtained produced important information from the theoretical point of view. 



In fig. |TB| the low-order diagrams iterated in the second RPA calculation are 
shown. The diagrams (a) are the usual RPA diagrams, while the other ones are 
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Figure 16: Representative low-order diagrams iterated in the second RPA theory. The 
(a) diagrams are the usual RPA diagrams. The (b),(c) and (d) diagrams are produced by 
the 2p-2h excitations. In addition to the self-energy insertions on a single particle line (b), 
there are ph linked diagrams of bubble (c) and ladder (d) type. 
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Figure 17: Longitudinal response at q=300 MeV/c in ^^C. The full line shows the full sec- 
ond RPA result. The dashed line has been obtained by neglecting the ph linked diagrams, 



the (c) and (d) diagrams of fig. 16 
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present only in a second RPA calculation. The (b) diagrams correspond to a cor- 
rection of the sp wave function, both hole and particle, produced by a self-energy 
insertion. The effects produced by these diagrams can be simulated by inserting an 
imaginary part in the MF potential. 

In fig. |l^ we present the ^^C longitudinal response at q=300 MeV/c obtained 
with a complete second RPA calculation (full line) and by neglecting the (c) and 
(d) diagrams of fig. |16| (dashed line). The similarity of the two curves in the QE 
region indicates that the most important diagrams beyond RPA are those of the 
self-energy insertion. This result is not valid in the giant resonance region | |Dro9U| ] 
as it is shown by the low energy part of fig. [l^. In ref. ||Dro89|| the effect of the (c) 



and (d) diagrams has also been investigated on the transverse response and it has 
been found to be slightly greater than in the longitudinal response. 

In spite of this difference the results of refs. |Pro87| , pro89|| indicate that the 
second RPA effects in the QE region can be simulated reasonably well by a complex 
optical potential. This approach has been used in MF | )Hor80| , phi89| , |Cap91| ] as well 
as in RPA |Pel85| , Bri87 , Bou89|, ^ag89 , Bou91 , Jes94 | calculations. The imaginary 
part of the optical potential removes part of the flux of the emitted particles and 
therefore this reduces the responses. The use of a complex, and energy dependent, 
optical potential has the disadvantage that the sp wave functions are no longer 
orthogonal, thus they do not form an orthonormal basis and cannot be used in 
many-body calculations like RPA. Furthermore the sum rules, controlling the total 
excitation strength, are no longer satisfled. These problems are usually neglected 
under the assumption that non-orthogonality effects are small. 

A different procedure to treat the FSI has been proposed in ref. ||Co'88|| . The 
basic assumption is to consider the self-energy of the system independent from the 
sp states: S(co') = A{uj) + iT{u)/2. The response containing FSI effects, R^^^ , can 
be obtained by a folding procedure: 

/•oo 

R^'^'iq, u)= dE R{q, E) [p{E, u) + p{E, -u)] (63) 

JO 

where the function p is related to the real and imaginary parts self-energy as follows: 



p{E,u) = — 



[UJ] 



(64) 



27r [E-u-/\{Lo)f + [V{uj)/2f' 

The r and A functions are flxed to reproduce the empirical values of the volume 
integrals of the optical potential | Mah82|| . In this treatment of the FSI the sum rules 
are conserved since the strength of the response is not reduced but it is redistributed 
on the energy spectrum. The general effect of FSI is a widening of the response width 
and a lowering of the peak value. 

In ref. [Po'88|| a q dependent nucleon effective mass m* such that m* jm < 1 , has 
been introduced to consider the non-locality of the static mean field. For a given 
m* /m value the following scaling relation holds: 



Rm* (g, ^) 



m" 



R 



m 



m 

g, — uj 

m 



(65) 
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Figure 18: ^^C responses obtained with CSM (dashed hnes). The full lines have been 
obtained applying the FSI corrections. 
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which shows that the response is lowered and that its niaxiniuni is shifted to higher 
energy. 

The inclusion of the effective mass still conserves the non-energy weighted sum 
rule, but it violates the energy weighted one. The size of this violation can be 
quantified by a relative enhancement factor which is of about 0.2 at g=300 MeV/c 
and 0.13 at g=550 MeV/c ||Co'88|| . 

The total effect of FSI and effective mass on the CSM responses of ^^C is shown 
in fig. [T^. These results are rather similar to those obtained using the optical 
potential. 

It is worth remarking that among the various corrections to the MF responses 
the FSI produces the biggest ones. It is curious, and unexpected, that effects re- 
lated to the many-particles many-hole excitations of the nucleus are more important 
than those related to Ip-lh excitations (like RPA). This fact can be understood by 
considering the strength of the residual interaction and the density of nuclear final 
states. The residual interaction in the QE region is very weak, therefore RPA ef- 
fects are small. On the other hand, the great number of 2p-2h states available in 
this region compensates the weakness of the residual interaction and this produces 
sizeable effects on the response. 



4 Comparison with experimental data 



Electron scattering experiments in the QE region have been performed at different 
laboratories for various nuclear targets. Longitudinal and transverse responses have 
been separated in ^H ||Ber871] , ^H [pi?87i PowSqi , ^He |[Ber87| , |Zgh94|| , ^^c [ |Bar83| 
^oCa piiS^i [MiJg^ , Piig^i piigBj lYit^Bj [Wiig7| , ^^Ca |Mez84 piiS3|, piiSq 
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Fe jMSOj |Mez84| , |Hot84| , |Mez85|| , ^ospb [|Zgh94|| and ^^^\J ||Bla8q . 



The comparison between theory and experiments is satisfactory in few-body 
systems, but for medium-heavy nuclei it is not possible to reproduce at the same 
time both responses. As example we show in figs. |19| and |2y the comparison between 
the results of our calculations with the available data on ^^C and ^°Ca. Each panel 
shows two curves, a dashed one corresponding to the standard treatment, i.e. CSM 
calculation with one body currents only, and a continuum curve containing all the 
corrections discussed in this paper excluding the relativistic effects. 

In the standard treatment the longitudinal responses are overestimated while 
the transverse ones are reasonably well reproduced. This fact provoked disconcern 
since the electromagnetic operator generating the longitudinal response, the charge 
operator, is supposed to be better understood than the current operator. 

A great theoretical effort has been done to explain the reduction of the longitu- 
dinal response. Various possibilities have been investigated, from relativistic effects 
[ Pel85| , [Bro894 |Weh89| , pior90|| to modifications of nucleon properties in the medium 
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Figure 19: Comparison of the ^^C experimental data of ref. [Bar83] with the CSM 
responses (dashed lines) and those obtained including the various corrections discussed in 
this paper, excluding the relativistic corrections. 
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Figure 20: The same as the previous figure for Ca. The data are from refs. |Mez84, 
Mez85|| (circles) and ||Yat93| , [Wil97|| (squared). 
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Nob81| , |Nob83| , Pro89b|| . It is now generally accepted that the FSI is the mecha- 
nism responsible for the lowering of the longitudinal response. The same effect acts 
also on the transverse response destroying the previous agreement (as is shown by 
the full lines of the figs. ^ and ^. The problem now has been moved from the 
longitudinal to the transverse response whose empirical values are larger than those 
predicted. 

A word of caution is necessary at this point. The responses are not directly 
measured but they are extracted from a set of cross section measurements. The 
extraction procedure is rather complicated and it is based upon various theoretical 
hypotheses on the structure and the behavior of the cross section. Some of the 
published data have been modified after a reanalysis [|Jou95| , |Jou96|| . The discussion 
on the reliability of the existing data and on the need of acquiring new ones is still 
open. It is not our intention to discuss this issue here. 



We would like to remark however that the new set of ^°Ca data [ |Yat93| , |Wil97| | 
(squares in fig. ^Oj) agrees rather well with our calculations. 



5 Conclusions 



In this article we have discussed the validity of the approximations made in what we 
have called the standard description of the QE excitation. This is a PWBA calcu- 
lation done within non relativistic MF nuclear model and using one-body currents 
only. 

The PWBA is usually adopted because it leads to a simple expression of the 
cross section, which allows the separation of longitudinal and transverse responses 
through a Rosenbluth procedure. Distortion effects produced by the Coulomb field 
of the nucleus are not very big. For nuclei with Z < 20 they can be rather well 
simulated by the effective momentum transfer prescription (|^). For heavier nuclei 
the focusing of the electron wave function onto the nucleus becomes relevant and 



DWBA calculations are necessary |Co'87a, Tra88 



A major source of uncertainty in the electron-nucleus interaction is related to 
the choice of the electromagnetic nucleon form factors |[Ama93a| . We have shown in 
fig. ^ the uncertainty band produced by the various possible choices. 

The size of the two-body currents effects is within this uncertainty band. We 
have studied two-body currents produced by the exchange of one pion. We have 
seen that the total effect of these MEC on the response is small because of the 
cancelation between the seagull and pionic current contributions [|Ama93b|| . Seagull 



and pionic currents modify the response in opposite directions and the final effect is 
rather small. We have also evaluated the contribution of the currents associated to 
the excitation of a virtual A which has been found to lower the transverse response 



Ama94a|| but to be rather small. 
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An estimate of the effects produced by relativistic corrections shows that their 
size is within the form factor uncertainty band as long as the momentum transfer 
value remains smaller than 500 MeV/c ||Ama96a , Ama96b , Ama98a . 

We started our investigation of the nuclear excitation mechanism by testing the 
need of using a finite system description for the QE excitation. We have compared 
CSM and FG responses and we found that they are rather similar if the value of the 
Fermi momentum is evaluated by considering the average nuclear density (see eq. 
(p7D). The agreement worsen by using the LDA ||Ama94b |. 



We have seen that RPA corrections are not important when a finite-range inter- 
action is used ||(Jo'88| , |Ama93b| . Zero-range interactions overestimate RPA effects, 
especially at high values of the momentum transfer. In infinite nuclear systems the 
choice of the residual interaction is inconsistent with the effective theory where it 
is used |[Bau984 [Bau98 b|| . From a pragmatic point of view these uncertainties are 
irrelevant since RPA effects are small. 

The largest correction to the standard description of the QE excitation is pro- 
duced by the FSI. Microscopic calculations done within the second RPA framework 
indicates that, in the QE region, the FSI can be rather well described in terms of an 
imaginary part of the optical potential [Pro87|| . We have developed an alternative 



approach to deal with FSI [|Co'88|] . Our results are similar to those obtained using 
the optical potential. The FSI enlarges the widths of the responses and lowers the 
peak values ||Ama93b[ . 

The comparison with the experimental ^^C data | )Bar83| | shows that, when the FSI 
corrections are included, the longitudinal responses are reasonably well reproduced 
while the transverse one are underestimated |IAma94al| . The same kind of calculation 



in ^''Ca does not reproduce the Saclay data ||Mez84|, [Mez85|| while the new MIT 



data [[Yat93| , [Wil971] are rather well described in both longitudinal and transverse 
responses. 

There are various problems still unsolved in the study of the QE excitation. At 
present the longitudinal responses seem rather well understood but the transverse 
ones are underestimates. There are exclusive data |[Ulm87| , |Ber90|| indicating that 
the two-nucleon emission in the QE region could be larger than that predicted by 
the MEG ||Ama93b[1 . This suggest the need of considering short-range correlations 
in the nuclear model ||Fan87t |Fab89| , [Lii90| , pH9l| , |Fab94| , |Fab97| , |Ama98b| , |Col 



Also the so-called dip region, the region between the QE peak and the peak of the 
A resonance, is at present not well described. 

The QE excitation of nuclei is a good testing ground of our knowledge of the 
electron-nucleus interaction and our understanding of the nuclear excitation mecha- 
nism. Since surface effects are not important |Ama94b| it is meaningful to compare 



directly the experimental data with the results of complicated many-body theories 
developed in infinite systems ||Fan8?l , |AIb89| , |Fab89| , |Cen94| , [Fab94| , |Uen97| , F^^b97 . 



We believe that the solution of the problems mentioned above and the interest 
of the theoreticians in comparing the results of their models with QE data require 
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new accurate and reliable data which could be rather straightforwardly obtained at 
the new electron accelerator facilities. 
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